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Holes in cubic semiconductors have effective spin 3/2 and very strong spin orbit interaction. Due 
to these factors properties of hole bound states are highly unusual. We consider a single hole bound 
by a spherically symmetric potential, this can be an acceptor or a spherically symmetric quantum 
dot. Linear response to an external magnetic field is characterized by the bound state Lande g-factor. 

We calculate analytically g-factors of all bound states. 
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I. INTRODUCTION 


Holes in cubic semiconductors originate from atomic 
p-orbitals. Therefore, the total angular momentum of a 
hole can take values j = 3/2 and j = 1/2. So the spin 
orbit split valence band consists of j = 3/2 and j = 1/2 
subbands. In semiconductors with large fine structure 
the j = 1/2 state is very high in energy and therefore is 
irrelevant. This is the case we consider here. Following 
tradition in the field, below we call the angular momen¬ 
tum of the hole “the effective spin” and denote it by 
letter S, j ^ S = 3/2. Studies of holes in cubic semicon¬ 
ductors have a long history. Pioneering theoretical work 
belongs to Luttinger who has applied the k ■ p method 
and the group theory analysis for hole quantum states 
near bottom of the band at the P pointi. Quadratic in 
momentum Luttinger Hamiltonian depends on Luttinger 
parameters. The Hamiltonian has a rotationally invari¬ 
ant part which is independent of the lattice orientation 
and it also contains a small correction which is related 
to the cubic lattice and which is not rotationally invari¬ 
ant. In the present work we neglect the small rotationally 
noninvariant correction and consider the rotationally in¬ 
variant Hamiltonian. 

The hole binding in an acceptor level or in a quan¬ 
tum dot is a problem of a tremendous experimental in¬ 
terest. Theoretically shallow acceptor levels have been 
studied thoroughly, from the spherical approximation^ 
to accounting the lattice corrections^, spin orbit split va¬ 
lence band^, and so-called central cell correctioit^^— . As 
the acceptor spectrum is quite sensitive to variations of 
Luttinger parameters, the experimental data have been 
used to determine the parameters^. The hole spectrum 
of the spherically symmetric quantum dot with infinite 
walls has also been calculated^. 

Zeeman splitting of acceptor levels in an external mag¬ 
netic field has been also studied theoretically. The most 
detailed numerical calculation includes the rotationally 
invariant approximation, cubic lattice correction, and 
also higher order magnetic field corrections^. Results of 
this calculation have been used in analysis of numerous 
experimental data, see also^^. Experimental data on the 
shallow acceptor spectrum in magnetic field have been 


obtained predominantly from infrared absorption spec¬ 
troscopy for different acceptors in GaAs^i^— , , 

Sji9-^ and other semiconductors. 

We already pointed out above that in the present work 
we disregard the spin-orbit split j = 1/2-band and limit 
our analysis by the rotationally invariant approximation 
in the Luttinger Hamiltonian. This implies that our anal¬ 
ysis is applicable for GaAs, InAs, InSb, and on the other 
hand the approximation is not justified in Si. We cal¬ 
culate analytically g-factors of hole bound states for ac¬ 
ceptors and spherically symmetric single hole quantum 
dots. At corresponding values of parameters our results 
are consistent with that of previous numerical computa¬ 
tions for acceptors^. 

Structure of the paper is the following. In Section H 
we discuss the Hamiltonian and approximations. In Sec¬ 
tion HI we calculate energy levels and wave functions 
of a hole bound by a spherically symmetric potential. 
This section is very similar to the paper—. Nevertheless 
we present this section because of two reasons, (i) it in¬ 
troduces the method and the notations which we use in 
following sections, (ii) We calculate previously unknown 
energy levels in a spherical parabolic quantum dot. Sec¬ 
tion IV presents our main result, analytical calculation of 
Zeeman splitting. In Section V we specially consider the 
“ultrarelativistic” limit, the limit when mass of the heavy 
hole is diverging. Section VI presents tables of g-factors, 
and our conclusions are summarized in Section VH. Alge¬ 
bra of tensor operators is described in Appendix A, and 
the derivation of spin-orbit radial equations is given in 
Appendix B. 


II. HAMILTONIAN 


Luttinger Hamiltonian for S' = 3/2 free hole in a zinc- 
blende semiconductor reads^ 

H,= {-n + ^ (vlSl + Plsl + plsl) 


(PxPy {Sa,, Sy } -f PyPz\^Sy, Sjg } “t" PzPx^^z , S^; } ) , 

mo 


X, y, z are the crystal axes of the cubic lattice, p is the 
quasi-momentum, m is the electron mass, 71, 72 and 73 
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are Luttinger parameters, and {...} denotes the anticom¬ 
mutator. This Hamiltonian can be rewritten as 


Hl= hi + ^72 


72 


CJ\2 




to Here we follow the standard conventions^ to use 
Hq from (21) with p —tt = p — eA. In this case the spin 
Zeeman interaction is 

Hz = -2K/isB • S, (4) 


where 


72 = 


272 + 373 
5 


where 2 k^bS is the spin magnetic moment, is Bohr 
magneton, and value of k depends on the materialSS. 
There is also cubic in spin Zeeman-like interaction 


The irreducible 4th rank tensor hjh depends on the ori¬ 
entation of the cubic lattice, the tensor is proportional to 
72 — 73 - Since in the large spin-orbit splitting materials 
72 ~ 73 the rotationally noninvariant part of the Hamil¬ 
tonian can be neglected. Hence the Luttinger Hamilto¬ 
nian can be approximated by the following rotationally 
invariant (independent of the lattice orientation) Hamil¬ 
tonian 


Hl^Ho = 


71 + 272 


5_ \ p 


5^-2^ (psr 

Zm m 


(1) 


Since different components of momentum commute be¬ 
tween themselves, Hq can be rewritten as 


Ha = 




( 2 ) 


n, = {S,,S,}--6,,-SiS+l) 
2 

Q ij {Pi; Pj } 2 ’ P 


H’z cx . (5) 

however it is small^^ and we disregard it. 

We restrict our analysis by linear in magnetic field ap¬ 
proximation. Therefore, gauge terms, p —>■ tt = p — eA, 
also result in a magnetic moment. As usually, see Refi^, 
we use the gauge A = [B x r]/2 to derive the orbital 
magnetic moment. All in all the operator of the total 
magnetic moment is where 

Mk = 2KSk + liLk- jijTijkTij. ( 6 ) 

Here the first term corresponds to @, the second term is 
the usual orbital magnetic moment which originates from 
the 7 i-term in Eq.(2|), and the last term is a spin depen¬ 
dent magnetic orbital moment which originates from the 
72 -term in Eq.([2]). The 3d-rank tensor in the last term is 

Bijk — 2 Pj} 4“ fo'fcn Pi}) : (7) 


The QijTij term in ([T]) represents the effective spin orbit 
interaction. Therefore, it is convenient to introduce the 
parameter 


and Tij is defined in ([ 2 ]). 

Thus, the total Hamiltonian which we consider in this 
work is 



71 


which characterizes the relative strength of the spin or¬ 
bit interaction. Eigenstates of the Hamiltonian o have 
definite helicity, the projection of spin on the direction 
of momentum, Sp = ±1/2 (light holes) and Sp = ±3/2 
(heavy holes). The corresponding masses are 

m m 

'rni = —yj——t = —vj-r . 3) 

7i(l±i/) 7i(l-i/) 

In the limit v ^ 1 mass of the heavy hole is diverging, 
nih —t oo, therefore we call this the “ultrarelativistic” 
limit, the limit corresponds to maximum possible spin 
orbit interaction. 

To account for a magnetic field B imposed on the sys¬ 
tem one has (i) to account for the vector potential via 
the long derivative p —>■ tt = p — eA, and (ii) to account 
for the Zeeman contribution oc S • B. Here e is the el¬ 
ementary electric charge. Since different components of 
7 r do not commute there is a well known ambiguity be¬ 
tween the gauge and the Zeeman terms. In particular, 
(tt-S)^ = |{ 7 rj, 7 rj}-{S'i,S'j} —|B-S. Therefore the Zee- 
man term added to Eq. (Ill) is different from that added 


H = Ho + U{r)-^lBm■^B) ( 8 ) 

where U{r) is spherically symmetric potential attracting 
the hole. 


III. ENERGY LEVELS AND WAVE 
FUNCTIONS OF A HOLE IN A SPHERICALLY 
SYMMETRIC POTENTIAL 

We already pointed out that this section is very sim¬ 
ilar to papers^i^. We present this section to introduce 
the method and the notations for further calculations of 
g-factors. Besides that we calculate previously unknown 
energy levels in a spherical quantum dot. The Hamilto¬ 
nian ([H]) is invariant under rotations around the center of 
symmetry of the potential. Therefore, while because of 
the spin orbit interaction the orbital angular momentum 
L = [r X p] is not conserved, the total angular momen¬ 
tum J = L ± S is conserved. The eigenstates of the 
Hamiltonian (l 8 |) are also the eigenstates of = J{J+ 1) 
and Jz = M. The Hamiltonian contains only the orbital 
scalar p^ and the orbital quadrupole Qij. Therefore an 
eigenstate can be a combination of at most two values of 
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orbital angular momentum with AL = 2, we denote by 
L the minimum of these two values 

= f{r)-\L,SJ,M)+9{r)-\L+2,SJ,M), (9) 

Radial wave functions are normalized 

OO 

f (I/p + |gp) = 1. If J > S' = 3/2, only two 
0 

values of L are possible: L = J — 3/2, L = J — 1/2. The 
special case of J = 1/2 realizes only at L = 1 or i = 2. 
In this case g(r) = 0 and the wave function is 


limit 1 / —>■ 1, see Eq.([3]). We will see that this results in a 
qualitatively different behavior of J > 3/2 and J = 1/2 
energy levels and g-factors in the limit ^ > 1 . 

To derive radial equations for functions / and g one has 
to use Wigner-Eckart theorem and algebra of spherical 
tensor operators. In Appendix A we remind the theorem 
and give definition of relevant spherical tensors. Equa¬ 
tions for / and g are derived in Appendix B, here we 
just present the equations. For J = 1I2 there is only one 
equation (we set h=l) 


d/(L, 5, J = 1/2, M) = f{r) ■ \L, S,J= 1/2, M), (10) 

Interestingly, while the states with J > S' = 3/2 are 
constructed of both heavy and light holes, the states with 
J = 1/2 are constructed of light holes only. Masses of 
light and heavy holes behave differently in the important 


1 

2mi 


L{L + 1) 


-f U{r) - E 


/ = 0 

( 11 ) 


We remind that here L can only take values L = 1,2. 


I 

For J > 3/2 and L = J — 3/2 the pair of coupled equations reads 


^ (l I 

2m V 2 J 


7i ( ^_v J-3/2 




/+ 


-f 


V3(J-l/2)(J + 3/2) 


2m 


2 J 


4m 
d, 


J 


. (j-i/2)(j + 3/2) 


2^^- (T+l/2)(J + 3/2) ,^^(^)-^ 


rr I '-'r 2 


1/71 v/3(J-1/2)(J + 3/2) 


drr — 


4m J 

Finally, for J >3/2 and L = J — 1/2 the pair of coupled equations is 


7i ( J + ^/2 


2m 


2 J+1 


dr 




-f 


t /71 V3(J-l/2)(J + 3/2) 


4m 


J + 1 


7i 

2 m 


1 + 


1 /J + 5/2 
2 J+1 


drr + —dr — 


(J + 3/2)(J + 5/2) 


1.71 V3(J-1/2)(J + 3/2) 


4m 


J+1 


drr - dr 


^2 

9+ 

(J + l/2)(J-3/2)' 


J.2 

/+ 


(J+l/2)(J + 5/2)' 


J.2 

-E 

9+ 

{J- 

- l/2)(J + 3/2)' 


9 = 0 


f = 0 


9 = 0 


f = 0 


( 12 ) 


( 13 ) 


Unlike Eq. (El) which contains only mass of the light hole, 
Eas. lfT^ and (fl^ cannot be expressed solely in terms of 
the heavy hole mass. Both heavy and light holes con¬ 
tribute in J > 3/2 states. However, it is obvious that in 
the limit —>■ 1 , m^ 00 , the heavy hole dominates in 
the wave function and hence the energy must scale as 


where a is the characteristic size of the wave function 
which also depends on v. Following Refj^ we will use the 
standard spectroscopic notations for the quantum states, 
for example l/S' 3/2 indicates that J = 3/2, the principal 
quantum number is n = 1 , and L = 0 (of course there is 
also L = 2 admixture). 

Numerical integration of Eqs. (in]),(ini),(ii3i) is straight¬ 
forward. Here we consider two special cases, (i) parabolic 
potential, and (ii) Coulomb potential. The parabolic po- 
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FIG. 1: Energy levels versus in a parabolic quantum dot. 
The levels with J = 1/2 are shown by solid red lines, 2Pi/2, 
3i3i/2, 4 Pi/ 2 from bottom to top. The levels with J = J> 12 are 
shown by solid black lines, IS3/2, 2P3/2, 35'3/2, 4P3/2, 3D3/2, 
4 F 3/2 from bottom to top. The levels with J = 5/2 are shown 
by dashed blue lines, 2 P 5 / 2 , 3 P 5 / 2 , 4 P 5 / 2 , 4 P 5/2 from bottom 
to top. The levels with J = 7/2 are shown by dashed-dotted 
green lines, 30^/2 and 4 P 7/2 from bottom to top. The level 
with J = 9/2, 4 P 9/2 is shown by the solid cyan line. Vertical 
dotted lines indicate values of v that correspond to particular 
compounds. 



FIG. 2: Radial wave functions of the ground state versus ra¬ 
dius. Black lines show /(r) and red lines show g{r). Solid lines 
correspond to GaAs and dashed lines correspond to InAs. 
The left panel correspond to the parabolic quantum well. The 
scale is ^ — l/yuum/^. The right panel corresponds to the 
Goulomb acceptor. The scale is the effective Bohr radius, 


tential we take in the form 

The ratio E/uj depends only on v. Plots of E/w ver¬ 
sus v for several lowest states are presented in Fig IT] 
The ground state is 153/2. While energies of states with 
J = 1/2 grow with v approximately linearly, the energies 
of states with J > 3/2 scale as E oc \/l — v in agree¬ 
ment with Eg. lfTTl) . For illustration we use Si, GaAs, 
and InAs. Values of Luttinger parameters from Ref^ 
are presented in Table |T| Our approach is valid when 
1 72 — 73 1 is small compared to 72 . Obviously, the ap¬ 
proach is not valid in Si, on the other hand it is quite 
reasonable in GaAs and even better in InAs. Values of 



71 

72 

73 

72 

V 

K 

Si 

4.285 

0.339 

1.446 

1.00 

0.47 

-0.42 

GaAs 

6.85 

2.1 

2.9 

2.58 

0.75 

1.2 

InAs 

20.4 

8.3 

9.1 

8.78 

0.86 

7.6 


TABLE I: Luttinger parameters for Si, GaAs, and InAs. 

1 / corresponding to Si, GaAs, InAs are shown in FiglTjby 
vertical dotted lines. In the left panel of Figl^jwe plot ra¬ 
dial wave functions of the quantum dot ground state for 
GaAs and InAs. The s-wave is dominating, the weight of 
the d-wave, 

00 

h = J drr‘^g‘^ (15) 

0 

is 21% in GaAs and 32% in InAs. 

Another example is attracting Goulomb center corre¬ 
sponding to acceptor impurity. 



where e is the dielectric constant. While this case has 
been thoroughly considered in the literature, we briefly 
present the results since we use them in the next Section. 
In this case the “atomic energy unit” is 

E* = . 

The ratio E/E* is plotted in FiglSj versus ly for several 
lowest states. In the right panel of Figl2] we plot the 



FIG. 3: Coulomb acceptor energy levels versus v. The level 
with J = 1/2, 2Pi/2, is shown by the solid red line, top panel. 
The levels with J = 3/2 are shown by solid black lines, IS 3 / 2 , 
in the bottom panel and 2 P 3 / 2 , 2 S 3 / 2 , in the top panel. The 
level with J = 5/2, 2 P 5 / 2 , is shown by the dashed blue line, 
top panel. Note that the vertical axis scales are different in the 
top and in the bottom panels. Vertical dotted lines indicate 
values of v that correspond to particular compounds. 

acceptor ground state radial wave functions for GaAs and 
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InAs. Again, the s-wave is dominating, the weight of the 
d-wave, Ii, is 26% in GaAs and 35% in InAs. Our results 
for the Coulomb acceptor are in excellent agreement with 

Rea. 


IV. LANDE G-FACTOR 


The operator of magnetic moment fOl is defined by 
Eq.®. Hence Lande g-factor of a quantum state with 
given J is defined by 

{JM\m^\JM) =g-M . 


The operator fOt contains a simple spin, oc S, and a simple 
orbital, oc L, contribution, as well as a more complex part 
dependent on the third rank tensor Ty ^, where the tensor 
is defined by (O. In order to apply the Wigner-Eckart 
formalism, we have to represent the third rank tensor in 
terms of irreducible tensors. It is easy to check that 


Tijk — g (^nikT^n % 

^ {^ik^j T ^jk^i ‘^^ij^k') 


where 


3 1 

T-f = ^ - -Sij{Xn,Pn} 


(16) 

(17) 


is the irreducible tensor of the second rank. For the mag¬ 
netic moment (0 we need only the convolution TijkTij. 
Therefore, we replace Tijk by the following simpler tensor 
which gives the same convolution. 


Tijk —t Tijk — r^^njkTiri ^jkLi- 


(18) 


Hence we can split the magnetic moment operator (|6]) 
into four parts = 9JIs -I- -I- 9112 + Slla which have 

different kinematic structures 


Ms = 2 kS 
niii = 7 iL 

V 

^2i = ll^TkiLk 

®l3i = -li^SnpTjfJrkj . 
The Lande factor splits correspondingly 
9 = Os + 9i + 92 + 93 


gs = 2 k 

91 = 7i 

92 = 7i 


(M 

M 

(L.) 

M 

V (n^Lj) 
4 M 


M 


93 = -7i- 


(19) 

( 20 ) 


where we average over the state \JM) given by Eq.® or 
Eg.lfTUl). 


A. 911s and Dili contribution 

The hrst two contributions, 9 s and gi, in Ea. (l20l) are 
usual textbook ones, they can be easily calculated via 
ordinary vector model^^: 

7 , S(S + l)-L(i + l)-4(L + |)/,'l 

-jpTi)-j 

7i f L(L -|- 1) — S(S -|- 1) -l- 4 (L -|- I) /i 

= -J(7TI)- 

Here Ii, Eg. dTSl) . is the weight of the component with 
higher orbital momentum L -|- 2. Notice that /i = 0 for 
J = 1/2. 



B. 9 JI 2 contribution 


In Eos. dTOll . d^ the operator 9112 is defined as product 
of Cartesian tensors, 9H2i oc 79^ = TikLk- On the other 
hand the tensor product theorem (IA6I) is formulated in 
terms of spherical tensors where the product is defined 
in terms of Clebsch-Cordon coefficients ic 


^9 Z^''^lq'2q-q'^q' ^q-q' 


(23) 


For Li and we use the conventions (IXTI) . (IA 2 I) to relate 
spherical and Cartesian components. Hence the relation 
between Cartesian and spherical components of D is de¬ 
termined by (l23l) . 


n(i) _ ^10 Tj 

~ *-^10,20-'^2: ~ y ■ 

Now we can use the general tensor product formula (IA6I) 
in combination with Wigner-Eckart theorem to find the 
matrix element of 


{L',S,J,M\Dz\L,S,J,M) = -M 


15 2J + 1 

T 7(7+1) 


X 


x \ L S J l(L||L(i)||L)-(5||r(2)||5)5iiq (24) 
[L S J j 


where {...} is the 9j-symbol and the reduced matrix ele¬ 
ments are 

(L||L«||L) = 7L(L + 1)(2L + 1) 

(5||t( 2)||5) = iv'(25- 1)2.5(2.5+1)(2.5 + 2)(25 + 3) 

O 

= 4\/5. (25) 

In Eg. (1241) we keep in mind that / and g-components of 
the wave function differ by AL = 2, and the operator Di 
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which is linear in L, does not connect these components. 
Hence, the g 2 -Lande factor takes the form: 

92 = 

= 71^ {{L, S, J, M\D,\L, S, J, M) • (1 - h)+ 

{L + 2, S, J, M\D,\L + 2, S, J, M) • h), (26) 

where the matrix elements are given by (HI- Substitut¬ 
ing 9j-symbols we arrive to the following answer 


J = 1/2 : 

biy 

92 = -7iy if L = 1 
92 =liv if L = 2 


J > 3/2 : 

V / J- 3/2 

z/ ( J + b/2 


2J+1 

J+l 

2J+1 

J 


h 

(1 - /i)) 


if L 
if L 


(27) 

J-3/2 
J - 1/2 


C. OTs contribution 


The calculation of the fOls contribution is very simi¬ 
lar to ‘0)12 ■ The operator OJls, Eos. (IT^ . (1201) . is defined 
as vector product of two irreducible Cartesian tensors r 
and of second rank, OJlsi oa Ri = SnijT^n Tkj ■ On 
the other hand the tensor product theorem (IA 6 I) is for¬ 
mulated in terms of spherical tensors where the product 
is defined in terms of Clebsch-Gordon coefficients 


The reduced matrix element of is given in (l25l) . Non¬ 
zero reduced matrix elements for tO) are 

X - 2 * 

+ = (31) 

X ^{r,Pr} + *(2T-l-1) 

X - i{2L+l)j , 

where Pr = —idr- The diagonal matrix element 
(L||T(2)||L) gives zero after the r-integration, so there 
are no // and gg terms in fDls. This is a direct conse¬ 
quence of the time reversal invariance, T^) changes sign 
under time reversal and on the other hand it is impossi¬ 
ble to construct a second rank tensor which is consistent 
with Wigner-Eckart and which changes sign under time 
reversal. This implies that g^ is zero for J =1/2 states. 
However, the fg contribution to SUla which originate from 
m is nonzero. A direct calculation with (1311) gives the 
following answer 

= (32) 

2(^L + 0/2+/3 , 

where 


= 7i^^ 


v'3(J-l/2)(J + 3/2) 


2J{J+l) 




_( 2 )y( 2 ) 

q' ^q-q' 


(28) 


For and we use the convention (IA2I) to relate 
spherical and Cartesian components. Hence the relation 
between Cartesian and spherical components of Ri is de¬ 
termined by (pS)) . 

RP = = ^R^- ( 29 ) 

According to (IA 6 I) . matrix elements of R^ break into 
the product of reduced matrix elements for and r: 


{L',S,J,M\rP\L,S,J,M) = 


J(J + 1) 


r 2 2 n 

xl S L' J WL'||r(2)||L). (5||r(2)||5'). (30) 

[ S' L J J 


/2 = J drr'^f{r)g{r) 

0 

GO 

h = j drr^ifg'- gf ) . (33) 

0 

We remind our convention, L is the minimum of two 
mixing orbital momenta, L corresponds to / and L -|- 2 
corresponds to g. 

Eqs. (l20l) . (1^ . (l22l) . (l27l) . and (1321) solve the problem 
of g-factor of the hole bound state. Answers for J = 1/2 
states are universal 

^’ 1 / 2 : 5 = yK - y(2-h 5i3) 

Di/2'■ 5 = —2zv-I- 71 ( 2 -I-z/) . (34) 

Answers for J > 3/2 states are not universal, they de¬ 
pend on radial wave functions. In Section VI we present 
the results for parabolic quantum dot. However, before 
calculating numerical values of the g-factor we discuss 
the “ultrarelativistic” limit z/ —>■ 1 . 
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V. G-FACTOR IN “ULTRARELATIVISTIC” 
u LIMIT 


The case of maximum possible spin-orbit coupling 
V ^ 1 (71 = 272 ) is important as a theoretical limit and 
also some real materials such as InSb and InAs are pretty 
close to this limit. As it has been mentioned above, en¬ 
ergies of J > 3/2 states in the limit ly —)■ 1 satisfy the 
approximate relation m- The states fall down on the 
bottom of the well. On the other hand, the states are 
described by equations dni) and (HSl). Because of col¬ 
lapse of the wave functions at = 1 we can neglect the 
energy and the potential in these Eqs. We keep only drr, 
^dr and ^ terms. Excluding the second derivatives by 
subtraction of two coupled radial equations, we find the 
following relation between functions f(r) and g(r): 

o^(f--f)= 9 ' + ^^^9 (35) 

\ r J r 

L = J-3/2, 

\ \l 3 (J+ 3 ^/ 2 j’ L = J— 1/2. 

Multiplying Ea. (l35l) by gr^ and separately by fr^ and 
then integrating by parts keeping in mind the normaliza¬ 
tion condition, f pr^dr = 1 — /i, one finds the following 
universal expression for the weight of the ^-component 
in the “ultrarelativistic” limit 


h 


3(J-l/2)/(4J), L = J-3/2, 
(J-l/2)/(4(J+I)), L = J-1/2. 


Interestingly, Eas. (l38)) are the same for L = J — 3/2 and 
L = J —1/2. Thus in the “ultrarelativistic” limit, ^ 1, 

g-factors of J > 3/2 states depend only on J and they 
are independent of a particular confinement (parabolic 
quantum dot. Coulomb acceptor,...) and independent of 
the radial quantum number. This universal behaviour is 
due to the collapse of heavy hole wave functions. Another 
point to note is the oc 1/[J(J -I- 1)] dependence of the g- 
factor (l38l) . This is similar to the g-factor dependence in 
diatomic molecules^^. Physical reasons for this are also 
very similar, in the molecules the dependence is due to 
the spin-axis interaction, for holes the dependence is due 
to (p • S)^ spin-momentum Luttinger interaction. 


VI. NUMERICAL VALUES OF THE LANDE 
FACTOR 


We have already pointed out that, in principle, Eqs. 
(Uni), (ED), ED, ED, dSD solve the problem of g-factor of 
the hole bound state. Still, to calculate numerical values 
of g-factors one needs to calculate integrals /i, / 2 , and 
I 3 . The integrals immediately follow from wave functions 
found in Section III. We have checked that our results for 
g-factors of Coulomb acceptor states perfectly agree with 
that found numerically and tabulated in Ref.—. Below we 
present our results for three lowest states in parabolic 
quantum dot, l/S' 3 / 2 , 2 P 3 / 2 , 2 P 5 / 2 , Tables [iTl HIH and 
Irrespectively. We tabulate 5 s / k , 51/71, 52/71, 33/71 
for different values of v. In the tables we also present 
separately the total orbital contribution g^ 


The 53 -contribution (1321) depends on the radial integrals 
I 2 and I 3 defined in Eg. (15^ . The relation (1551) is not 
sufficient to determine I 2 and Is separately. However, 
remarkably, the combination 2{L -|- 3 / 2)/2 -|- Iz which en¬ 
ter Eq. ED can be found from (ED. Using the same 
derivation as (l36l) we find 


2{L + 3/2)/2 + h 



Substitution of the integrals (IMl) and (1^ in Eqs. (l20l) . 
(ED, (ED, (ED, a-nd (IMll gives the following values of the 
g-factors for J > 3/2 


9k 


2J(J-|- I) 


9s = 

91 = 7i (^1 - 

7i 

92 = 

93 = 

9 = 


4J(J-f-I) 
9 


4 4J(J-hI) 


1 - 


371 _ 

4 4J(J-fI) 

9(2k - 7i) 

4J(J-hI) ■ 


( 38 ) 


5l = 3i+32 + 53, g = 9s+gL- 

All lines in each table except of the lowest line are calcu¬ 
lated with Eqs. ED, ED, ED, ED, ED and with in¬ 
tegrals /i, / 2 , I3 calculated with functions / and 5 found 
numerically in Section III. The lowest lines are calculated 
directly with the “ultrarelativistic” Eqs. (ED- 


V 

9s/>=i 

51/71 

52/71 

53/71 

5n/7i 

0.6 

1.83 

0.087 

-0.052 

-0.285 

-0.25 

0.65 

1.78 

0.109 

-0.071 

-0.337 

-0.299 

0.7 

1.73 

0.135 

-0.095 

-0.392 

-0.352 

0.75 

1.67 

0.166 

-0.125 

-0.449 

-0.408 

0.8 

1.59 

0.203 

-0.162 

-0.505 

-0.465 

0.85 

1.51 

0.246 

-0.209 

-0.556 

-0.519 

0.9 

1.408 

0.296 

-0.266 

-0.596 

-0.566 

0.95 

1.3 

0.349 

-0.332 

-0.615 

-0.597 

1.0 

1.2 

0.4 

-0.4 

-0.592 

-0.6 


TABLE II: g-factor contributions for IS 3/2 ground state in 
the parabolic quantum well at different values of v. 


In all the cases there is a very strong compensa¬ 
tion between spin, gs, and orbital, 5 l, contributions 
to the g-factor. For example, in GaAs (zz = 0.75), 
gisy, = 2.0 - 2.8 = - 0 . 8 , g2P,^, = 1.6 - 2.6 = - 1 . 0 , 
32 P 5/2 = 0-37 — 0.40 = -1-0.47. These compensations are 




















V 

5s/k 

51/71 

52/71 

53/71 

5i/7i 

0.6 

1.375 

0.312 

-0.293 

-0.251 

-0.231 

0.65 

1.362 

0.319 

-0.313 

-0.288 

-0.282 

0.7 

1.347 

0.326 

-0.332 

-0.327 

-0.333 

0.75 

1.331 

0.334 

-0.349 

-0.369 

-0.384 

0.8 

1.314 

0.343 

-0.365 

-0.412 

-0.434 

0.85 

1.294 

0.353 

-0.380 

-0.457 

-0.484 

0.9 

1.271 

0.365 

-0.392 

-0.504 

-0.532 

0.95 

1.243 

0.379 

-0.400 

-0.553 

-0.574 

1.0 

1.2 

0.4 

-0.4 

-0.592 

-0.6 


TABLE III: g-factor contributions for the first excited 2 P 3/2 
state in the parabolic quantum well at different values of v. 


V 

5<,/k 

51/71 

52/71 

53/71 

5i/7i 

0.6 

0.897 

0.552 

-0.016 

-0.427 

0.108 

0.65 

0.841 

0.579 

-0.045 

-0.483 

0.052 

0.7 

0.784 

0.608 

-0.078 

-0.534 

-0.004 

0.75 

0.728 

0.636 

-0.116 

-0.578 

-0.058 

0.8 

0.674 

0.663 

-0.155 

-0.615 

-0.108 

0.85 

0.625 

0.687 

-0.197 

-0.644 

-0.153 

0.9 

0.582 

0.709 

-0.237 

-0.665 

-0.193 

0.95 

0.545 

0.728 

-0.277 

-0.679 

-0.228 

1.0 

0.514 

0.743 

-0.314 

-0.677 

-0.257 


TABLE IV: g-factor contributions for the second excited 
2 P ^/2 state in the parabolic quantum well at different values 
of V. 

somewhat similar to the compensation found in a partic¬ 
ular situation in Refi^. Here we see that the compen¬ 
sation is a generic effect related to the “ultrarelativistic” 
formula (|38)l . 

VII. CONCLUSIONS 

We have developed analytical tensor theory of Lande 
g-factors in spherically symmetric bound states of holes 
in cubic semiconductors. The Lande factors are given by 
Eqs. (Uni), (EH), (EH, (EH, (EH integrals Ii, / 2 , I 3 

which depend on the quantum state wave function, see 
Eqs. El), EH, (EH- We have shown that g-factors only 
weakly depend on the type of confinement (Coulomb 
acceptor, parabolic quantum dot, ....). Behaviour of 
g-factors in the “ultrarelativistic” limit is universal and 
completely independent of the type of confinement. 
This “ultrarelativistic” behaviour enforces a strong 
compensation between spin and orbital contributions 
to g-factors. It is worth noting that the developed 
tensor technique provides a solid basis for calculation of 
g-factor of asymmetric quantum dots, this problem will 
be addresses separately. 
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Appendix A: irreducible tensors, Wigner-Eckart and 
tensor product theorems 


culations. All tensors are defined in three-dimensional 
space. A tensor of the rank i is irreducible if 

its components can be transformed over some irreducible 
matrix representation of the rotational group SO(3). 
Consequently, the number of independent components of 
Riii 2 ...it coincides with the dimension of the representa¬ 
tion dimS)f = 2 £+ 1, whence we may choose the basis of 
these components Rm , rn = —£,... ,£ which have been 
transformed exactly like spherical harmonics This 
basis is usually called spherical basis. Spherical compo¬ 
nents of the vector R can be related to its Cartesian 
components in the following way 

rW ^ ± . (Al) 

Similarly spherical harmonics of the second rank irre¬ 
ducible tensor can be related to its Cartesian com¬ 
ponents as 

Rq^^ = Rzz, R±l = T (Rxz ± iRyz) , 

-^±2 = - Ryy ± ‘^iRxy) , (A2) 

The Wigner-Eckart theorem reduces the matrix ele¬ 
ment of a spherical tensor between states with given an¬ 
gular momentum \J,Jz = M) to the 3j-symbol2^ 

(A3) 

Importantly, the reduced matrix element {J'\\R^^'>\\J) is 
independent of the projections M and M'. 

Consider tensor product of two irreducible tensors 
and Q*)™) of ranks k and m respectively. The tensor 
product is also an irreducible tensor of rank \k — m\ < 
( < k + m. It is defined as 

R(o = (r('=) X ^ (A4) 

q' 

where C is the Clebsch-Gordon coefficient. The scalar 
product is a special case of the tensor product 

k 

(j’lfc) . gW) = (A5) 

q— — k 

If two irreducible tensors T and Q in the tensor product 
are independent, i.e. they act in different sub-spaces (say 
spin and orbit) then the matrix element of the tensor 
product is disentangled via the product of reduced matrix 
elements of T and Q: 


In this Appendix we summarize the major equations 
from the SO(3) tensor algebra^^ which we use in our cal- 
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(j( j'j'M'|(r« • g«)| Ji J 2 JM) = 5j,j ■ 5 m'm{-i) 


J+J1+J2 j j[ Ji k 

J2 J2 J 


fc m I 

J[ J!, J' i(J{||rW||Ji)(J'||QW||J2)(A6) 

Ji J2 J J 

(j(||rW||Ji)(j'||gW||j2), 


(A7) 


r 


where = y/{2( + 1)(2J + 1)(2J' + 1), and curly 

braces denote 9j and 6j symbols. 


Appendix B: spin-orbit radial equations 

To derive equations for radial functions of the bound 
hole we have to calculate matrix elements from tensors 
T^^'> and defined in ©• Since the general dependence 
of projections is given by the Wigner-Eckart theorem, Eq. 
(IA31) . we need to know only reduced matrix elements. 
Let us calculate for example the diagonal matrix element 
(L||g(^^||L) using quantum states \L,Lz) = \L,L) with 
Lz = L. According to Wigner-Eckart 




L 2 L 
-L 0 L 


(Bl) 


where = Qzz =‘^{pl — and the 3j-symbol i 


iO 

L 2 L 
-L 0 L 


IS 


L{2L- 1) 


(L + l){2L+l){2L + 3)' 


(‘2') 

To calculate (L,L|Qq |L, L) we use the explicit form of 
p1 in spherical coordinates: 

p1 = — f cos^ 9 ■ drr + -dr ■ (sin^ 9 — sin(20) • dg)+ 


—j(sin^ 9 ■ dee + sin(26») • de) 


The spherical harmonic = \L,L) has the following 
well-known form 




Performing the straightforward integration over the an¬ 
gles we find the matrix element 


(L,L|Qf |L,L) = 


4L 


3(2L-H3) 




and hence using (ED we find the reduced matrix ele¬ 
ment. The procedure to find the offdiagonal matrix ele¬ 
ment (L-|-2||g*^^^||L) is absolutely similar. All in all this 
gives 




(B2) 


X ( 5 ,, + -dr - 


LiL + 2) 


a 2 L + 1 

(Jrr - dr 


(L-2||<3»||L) = - 2 ^|L^ 


X i drr + 


2L + 1 


L^-1 


The reduced matrix element of is presented in 
Ea. (l25l) . Substituting the wave function in the 
Schrodinger equation with Hamiltonian ([5]) and taking 
the projections on spherical harmonics \L, S, J, M) we 
finally obtain equations (HD), (na), and (na for radial 
functions. When calculating matrix elements of QijTij 
in Eq. ([a one has to remember the following relation be¬ 
tween the Cartesian dot product and the spherical dot 
product (IAtI) . QijTij = |(Q • r). 
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